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Abstract 

Context. The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transi¬ 
tion disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. 
Aims. We aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and 
understand their origin. 

Methods. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 pm 
to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. 
Results. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non- 
axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO 
are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the 
star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model 
requires a high disk aspect ratio (H/r ~0.20 at the planet locations) to account for the large pitch angles which implies a very warm 
disk. 

Conclusions. Our observations reveal the complex morphology of the disk MWC 758. To understand the origin of the detected 
features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed. 

Key words. 


1. Introduction 

The planetary systems detected so far show a great diversity in 
their architectures and natures which may originate in the variety 
of initial conditions of their formation ( [Mordasini et al.|2012| . 
It is thus of fundamental importance to study the physical con¬ 
ditions in the circumstellar disks in which they form, and the 
processes that rule their evolution. Direct-imaging observations 

* Based on observations performed with SPHERE/VLT under pro¬ 
gram ID 60-9389(A) 


at high angular resolution of young disks provide an efficient 
way to constrain these processes because they may leave clear 
signatures in the disk structure. As an example, planet-disk in¬ 
teractions will naturally generate asymmetric features, such as 
vortices, spiral arms or eccentric gaps, and may be used to in¬ 


directly detect planets embedded in the ir parent disk (e.g Kley 
|& Nelson]|2012[ |Baruteau et al.|[20T4] ). In this context, disks 
with gaps or cavities, called transition disks, are prime targets 
to study because they may be at an advanced stage of disk evo¬ 
lution ( Espaillat et al.|2014| . Asymmetric features have already 
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been observed in transition disks with high-contrast observations 
in the infrared (IR), which probe the disk surface layers (e.g., 


Hashimoto et al. 2011; Muto et al. 

2012; 

Rameau et al. 2012 

Casassus et al. 2012[[Boccaletti et a 

12013 

JCanovas et al. 2013 

Garufi et al. 2013:; Avenhaus et al. 2014; T 

halmann et al. 2014). 


Some of these features wer e traced deeper in the di sk with CO 
submillimeter imaging (e.g., Christiaens et al .||2014| ). 


In this context, the transition disk around the Herbig A5 
star MWC 758 (HD 36 112) i s particularly interesting. It is 
3.5+2 Myr old (Meeus et al .|2012 ), and is located near the edge 


of the Taurus star-forming region. The revi sed Hipparcos paral 
lax data provide a distance of 279+^j pc ( van Leeuwen||2007 ) 


It is most likely a single star, as any stellar companion down 
to a mass limit of 80 Mj up was ruled out based on sparse aper¬ 
ture masking observations ( [Grady et al.||2013| . Submillimeter 
observations show a circumstel lar disk extendin g up to 385+26 
astronomical units (au) ( [Chapillon et al.||2008| ), with an incli¬ 
nation and position angle of 21° +2° and 65°+7°, r espectively, 
as determined from the CO line ( [Isella et aLj[2010 ). Modeling 


indicates that the surface density steeply increases between 40 
and 100 au, which supports the idea that the disk might be grav¬ 
itationally perturbed by one or more companions. The observa¬ 
tions marginally resolved a cavity with a radial extent in the 
subm i llimeter estimated be tween 70 and 100 au ( [Isella et~aL| 
20 10[ [Andrews et aL||2011| ). Using Ks (2.15//m) direct imag¬ 


ing and H -band (1.65 yum) polarimetric imaging with the High 
Contrast Instrument for the Subaru Next Generation Adaptive 
Optics (HiCIAO), [Grady et"aL| ( |2013[ ) found no evidence for a 
cavity in scattered light down to 0.1" (28 au) from the star and 
detected two small-scale spiral arms. The authors suggested that 
the southeast spiral might be launched by a 5+ 2 M Jup companion 
at 1.55" (432 au; beyond the submillimeter dust emission) in a 
dynamically warm disk. 

In this Letter, we report IR polarized intensity observations 
of the transition disk MWC 758 obtained with the SPHERE in¬ 
strument in the Y band (1.04/im).The use of an extreme adaptive 
optics system, combined with high-contrast coronagraphic and 
differential imaging techniques led to deeper observations, at a 
smaller inner working angle (IWA; 0.093", 26 au) and a higher 
angular resolution (0.027", 7.5 au) than previously achieved. 


tribution of small dust grains in the disk surface layers through 
their scattered light by suppressing the unpolarized stellar light 
(e.g. |Apai et aL||2004| |Quanz et aL||2011| ). IRDIS provides two 
beams, in which wire-grid polarizers are inserted in DPI mode, 
and lead to ordinary and extraordinary polarization states. The 
half-wave plate (HWP) that controls the orientation of the po¬ 
larization was set to four positions shifted by 22.5° in order 
to construct a set of linear Stokes vectors. To minimize uncer¬ 
tainties related to the relative positioning of the star with re¬ 
spect to the coronagraphic mask, the frame registration was ac¬ 
curately performed using a background star located at -2.35" 
from MWC 758. The data were corrected for the distortion and 
true North. We then followed a procedure developed for high- 
contrast imaging and used the double-ratio method to calculate 
the Stokes parameters Q and U (Avenhaus et al. 2014| . Since the 
scattered light from a circumstellar disk is expected to be linearly 
polarized in the azimuthal direction under the assumption of sin¬ 
gle scattering, it is beneficial to describe the polarization vector 
field in polar rather than Cartesian coordinates. We define the 
polar-coordinate Stokes parameters Q#, as: 

Qcf> = +<2cos(2®) + U sin(20) (1) 

= -Q sin(2d>) + U cos(20), (2) 

where (f> refers to the azimuth in polar coordinate, <D is the po¬ 
sition angle of the location of interest (x, y) with respect to the 
star location (xo,yo), and is written as: 


X- X() 

® = arctan-+ 6. 

y-yo 


(3) 


6 corrects for instrumental effects such as the angular mis¬ 
alignment of the HWP. In this coordinate system, the az- 
imuthally polarized flux from a circumstellar disk appears as 
a consistently positive signal in the Q$ image, whereas the U<p 
image remains free of disk signal and provides a convenient es¬ 
timate of the residual noise in the Q$ image ((Schmid et al.|2006J). 
Because the DPI mode of IRDIS was still in commis"sioning at 
the time of publication, the absolute level of polarization could 
not be accurately retrieved, and we are only able to report polar¬ 
ized intensity images with arbitrary surface brightness levels. 


2. Observations and data reduction 


The observations were carried out on 2014 December 5 and 2015 
March 3 with the SPHERE instrument ( [Beuzit et al.|2008[ Beuzit 
et al. in prep.), equipped with an extreme adaptive-optics (AO) 


system ( [Fusco et al.|2006[|Petit et al.|2014|[Sauvage et al.|2014 ) 
at the Very Large Telescope at Cerro Paranal, Chile. The obser¬ 
vations were executed in Science Verification Time. MWC 758 
was observed in the F-band filter (1.04 /mi) using the infrared 
dual-band imager and spectrograph (IRDIS; Dohl en et al.|2 008; 
Langlois et al.|2014) , with a plate scale of 12.25 milli-arcseconds 
(mas) per pixel. We used a 185 mas-diameter coronagraphic 
mask combined with a Lyot apodizer. MWC 758 was observed 
during -51 minutes at both epochs. These data were taken in 
challenging AO conditions with moderate to poor seeing (be¬ 
tween 0.9" and 1.2" in the optical), high airmass (between 1.5 
and 2), and relatively short wavelength. The Strehl ratio reached 
about -70% in the H band, that is, -45% in the Y band. The 
December dataset suffered from a number of defects (regard¬ 
ing spider alignment and centering onto the mask), leading to a 
lower signal-to-noise ratio (SNR) than in the March dataset. 

The images were obtained using differential polarimetric 
imaging (DPI), a powerful technique for studying the spatial dis¬ 


3. Polarized intensity images 

Figurefllshows the Q$ images at both epochs and the radial map¬ 
ping ofthe March data, obtained after deprojection (with i = 21°, 
PA = 65°;[isella et al. 2010). Each pixel was scaled with r 2 to 
compensate for the r~ 2 dependency of the stellar illumination. 
The dark central region corresponds to the area masked by the 
coronagraph. The residual spiders from the coronagraph are vis¬ 
ible in the December 2014 image (Fig. |T] left). Scattered light is 
detected from -0.093" to -0.53" (-26 to -148 au) in polarized 
intensity, and we did not resolve any depleted cavity. Six distinct 
non-axisymmetric features are detected, the innermost ones be¬ 
ing at previously inaccessible radii. The features are indicated 
by arrows in Fig. [T] (1) a spiral arm located in the south and 
east (hereafter, the SE spiral) at deprojected distances r -0.26"- 
0.53"; (2) an arc in the west at r ~0.2"-0.3"; (3) a spiral arm 
in the northwest (hereafter, the NW spiral) at r ~0.22"-0.34"; 

(4) an arc located north of the NW spiral at r ~0.32"-0.53"; 

(5) substructures inside the spirals at r ~0.1"-0.2" and (6) a sub¬ 
structure branching off the SE spiral at r ~0.21"-0.27". We note 
that although there is some clear scattered light signal very close 
to the coronagraph, the exact morphology of feature (5) may be 
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Figure 1. Left and middle: polarized intensity images (Q^) obtained in December 2014 and March 2015, respectively. East is toward 
the left. Right: radial map of the deprojected Q$ image from March 2015 using i = 21° and PA = 65°. The dashed line indicates a 
radius of 0.23". Each pixel has been scaled with the square of its distance from the star, r 2 , to compensate for the r~ 2 dependency 
of the stellar illumination. The color scale is arbitrary. 


affected by the efficiency of the AO correction and the corona- 
graph centering. App endix |A| provides the Uq images (Fig.lA.l} 
and SNR maps (Fig. |A.2| ). Tne features are detected above the 
disk background at more than 3-<x and do not appear in the U<p 
images. This suggests that they are real and not diffraction resid¬ 
uals that would affect the Q$ and images alike. Although 
one could naturally see a single spiral arm in the west, extending 
from PA ~ 180° to 350°, the radial mapping (Fig. [I] right) shows 
that it consists of different features that have a very distinct de¬ 
pendence on the PA , with a sharp transition at PA ~ 240°. In 
the December 2014 image, feature (2) hosts an area that is ~ 1.5 
times brighter than the NW spiral, at PA ~ 240° - 260°. In the 
March dataset, the area at PA ~ 260° in contrast seems fainter 
than the background spiral. Since it is very close to the resid¬ 
ual spiders, it is difficult to conclude whether it is an instrumen¬ 
tal effect, or if it has a physical origin which, considering the 
timescale of the variations, would point to a variability in the in¬ 
ner disk. While the detected structures agree with the HiCIAO 
//-band image ( |Grady et al.||2013 ), we did not detect a signifi¬ 
cant difference in the extent of the scattered light signal between 
the west and the east sides. 


O(r) 


O r + 


sgn(r - r c ) ( r 

h r 


i +P 


\+P 1 


~—(- 
a+P\r c 


sgn(r -r c )( 1 




\+P 1-a+p 


(4) 


4. Spiral feature modeling 

Planets embedded in circumstellar disks are known to launch 
spiral waves at Lindblad resonances both inside and outside of 
their orbit (e.g. Ogilvie & Lubow|2002| ). The shape of the spiral 
wake is determined by the disk rotation profile and by the sound 
speed (i.e., temperature) distribution in the disk. In this section, 
we attempt to fit the shape of the two regions where there is a 
significant departure from circular symmetry, features (1) and 
(3), with a model based on the spiral density wave theory. As 
the spiral features are very similar in both datasets, and since the 
observations in March were obtained under better conditions, we 
only fit the corresponding dataset. We assume that the observed 
scattered light traces small dust grains that are well coupled to 
the gas, and thus, that it indirectly traces the gas. Following the 
prescription of |Rafikov| ( |2QQ2| ) and |Muto et al. ( |2012| ), the spiral 
wake in polar coordinates (r, <D) follows 


where a and p are the exponents of the disk rotation and 
sound speed profiles, respectively: D oc r ~ a , c s oc r~h . h c = 
c s (r c )/r c £l(r c ) is the ratio of the pressure scale height H to the 
radius (also called disk aspect ratio) at the location of the planet 
(r c , O c ). The pitch angle of the spirals depends on the disk tem¬ 
perature (hence on the aspect ratio) and the distance from the 
launching planet. The flaring index, a - ft - 1, determines the 
variation of H/r with radius. Equation [4] is valid in the linear 
or weakly nonlinear regimes, in which a single spiral wave is 
launched by an embedded planet, and approximates well the 
shape of the density wave given by the WKB theory ( Rafikov 

|20fe| . 

We attempted to fit both spirals simultaneously assuming 
that they are launched by two planets in circular orbit at different 
radii in the disk. We assumed that the disk is in Keplerian rota¬ 
tion and fixed a= 3/2. Varying ft has little influence on the fit, and 
we fixed /3=0.45, following Andrews et al. ( 201 1| ). We consid¬ 
ered that h c is a global disk parameter, meaning that the values 
of aspect ratios at the locations of the planets should be con¬ 
sistent. This leads to five free parameters in total. We restricted 
our models to planets located inside the submillimeter cavity 
whose extent is -0.36" wi thin large uncertainties ( Isella et al. 
2010; Andre ws~et al.|20lf| and considered disk aspect ratios of 


at most 0.20 to be consistent with the modeling of the spectral 
energy distribution (SED) ( [Andrews et al.|20lT} . The parameters 
are thus varied as 0.03 < h c < 0.20, 0.01" < r C5N w/sE < 0.36", 
0 < 0 c ,nw/se < 360°, in 20, 50, and 50 linearly spaced values, 
respectively. 

We deprojected the image and fit the locations of the surface 
brightness maxima along a set of azimuth angles (PA ~70°- 
200° (80 linearly spaced values) and 280°-340° (50 values), 
for the SE and NW spirals, respectively). We minimized a x 1 
function assuming that the uncertainty on each location is the 
FWHM of the point spread function (PSF). The best-fit param¬ 
eters are for the SE arm: (h c , r c , O c ) = (0.20, 0.253", 72°), and 
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slightly higher than the predictions from SED and submillimeter 
continuum modeling (-200K at these locations; Andrews et al. 


2011), but significantly different from the temperature measured 
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with CO (~53 K at 50 AU; |Isella et al.|2010| >. Such a high disk 
aspect ratio also reduces the contrast and the observability of the 
spirals in scattered light ( [Juhasz et al.|2014| ). This suggests that 
the vertically isothermal assumption implied in the spiral model 
may not hold (e.g. |Richert et al.|2015| . Alternatively, a planet on 
an eccentric orbit, an inclined planet or a planet massive enough 
to exert disk eccentricity may affect the shape of the spirals and 
induce a different pitch angle than in the case of a circular orbit. 

The two spiral arms were interpreted in the context of planet- 
disk interaction, but other mechanisms, such as self-gravity, 
can trigger spiral features with lower winding numbers (Rice 


et al. [2006b] ), while keeping enough contrast to detect them in 


Figure 2. Best-fit model in dashed lines. The locations of the 
planets are indicated by white diamonds. The dotted ellipse is a 
projected circle with radius 0.25". The color scale is arbitrary. 


for the NW arm: (h c , r c , O c ) = (0.20, 0.196", 266°). This cor¬ 
responds to planets located at -55 and 71 au for the NW and 
SE arms, respectively. The spiral shapes with these parameters 
are plotted in the polarized intensity image in Fig. [2] As noted 
by Muto et al. ( 2012] ), the degeneracy of the model is signifi¬ 
cant. Figure [BT shows the 2D probability distributions for each 
pair of free parameters. The Bayesian probabilities were derived 
as exp(-^ 2 /2) and normalized such that their integral over the 
whole parameter domain is unity. We obtained a lower limit of 
the disk aspect ratio h c = 0.2. A high value of h c is needed to ac¬ 
count for the large pitch angle, and while higher values may be 
favored by the spiral fit, they would produce too much infrared 
flux to account for the SED. In addition, we note that h c and r c 
are not independent. The lower h c , the closer to the spiral the 
planet must be, as the spiral pitch angle is only large very close 
to the planet. 

Although we restricted the planets to be inside the submil¬ 
limeter cavity, we note that the only way to fit the northern tip of 
the SE spiral (up to PA -45°) is to consider a planet outside of 
the submillimeter continuum emission as in |Grady et al. ( 2013 ). 
However, extending the range in r c to 2" strongly increases the 
degeneracy of the model. We do not favor this solution because 
companions are likely to be located inside the submillimeter cav¬ 
ity since other clearing mechanisms, such as photo evaporation, 
can be ruled out by near-IR interferometric observations of the 
inner disk ( [Isella et al.|2008[ ). 


scattered light (Pohl et al, in prep.) and at longer wavelengths 
(Dipierro et al.| 2014] ) . The disk mass has been estimated to 
be -0.01 M 0 from submillimeter observations (Andr ews et al.| 
201 1| ), which is probably too low to trigger gravitational insta¬ 


bilities, although these estimates are very uncertain and strongly 
depend on the assumed dust opacities and gas-to-dust ratio. If a 
significant part of the solid mass is in planets or particles larger 
than 1 cm, the gas mass can be much higher than currently esti¬ 
mated. 

Assuming that the observed spirals are trailing, the SE spi¬ 
ral arm is located on the far side of the disk and interestingly, 
does not appear darker than the near side (unlike SAO206462, 
Garufi et al. 2013| ). This may indicate an isotropic scattering 


5. Discussion and conclusions 

We presented 7-band polarized intensity images of the scattered 
light from the protoplanetary disk around the Herbig Ae star 
MWC 758 and studied the morphology of non-axisymmetric 
features. The image shows six distinct features, including arc¬ 
like features and two spiral arms with a large winding angle. 
Fitting the shape of the spirals, with a model based on the spiral 
density wave theory, showed that only a warm disk with a high 
aspect ratio (~0.2) at the locations of the planets can lead to such 
large pitch angle spirals ( |Rafikov|2002| . Assuming that the disk 
is in hydrostatic equilibrium, this corresponds to midplane tem¬ 
peratures of -315 K and -250 K at 55 and 71 au, respectively, 


phase function of the dust grains, which would imply submicron 
sized grains, or a high polarization efficiency. We note that the 
detected scattered light is quite symmetric and the western side 
of the image does not show the extended emission seen in the 
HiCIAO image ( [Grady et al.|2013| ), which was interpreted as be¬ 
ing due to an asymmetric irradiation of the disk. Considering that 
the observations were obtained only three years apart, it is un¬ 
likely that this would result from variable shadowing and should 
be investigated further. 

We found no indication of a fully depleted cavity in micron- 
size dust grains beyond the coronagraph radius (consistent with 
12 CO peaking at the star position; Isella et al. 2010 ) while a cav¬ 
ity in millimeter grains has been marginally resolved at a much 
larger radius (Isel la et al.|2010[ [Andrews et al.|201l] ). Such spa¬ 
tial segregation of small and large grains can be a natural out¬ 
come of particle trapping at the edge of a cavity carved by a 
planet (Pinilla et al. 2012, 2015]). The requirement of a contin¬ 
uous replenishment of small particles through the cavity trans¬ 
lates into a maximum planet mass, above which any companion 
would filter all dust particles ( [Rice et alT 2006a[|Zhu et al.|2012 ). 
With typical disk viscosities, a companion with a planet-to-star 
mass ratio above 10 -2 (i.e., a 5.5 Mj up planet around MWC 758) 
would filter all dust grains ( [Pinilla et al.|2012| ). 

Features (2) and (3) have a very distinct dependence on the 
PA , with a sharp transition at PA -240°, therefore no spiral 
model can fit the full NW structure (features (2) and (3) to¬ 
gether). A projected circle of 0.25" radius is shown in Fig. [2] 
and can only partly account for the region between PA -190- 
270°. The circular region between features (2) and (3) may trace 
the edge of a cavity, although high-resolution submillimeter ob¬ 
servations are required to determine whether this is the case. 
Alternatively, a slightly eccentric ellipse (e~0.1) oriented with 
PA=- 30°, that would trace an eccentric gap, can fit part of fea¬ 
ture (3) in addition to feature (2). In addition to the two spiral 
arms, four more structures are detected. Feature (4) is located 
beyond the NW spiral and indicates that the NW arm does not 
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shadow the outer disk in this direction. Its small angular extent 
is puzzling, and its nature is unclear. Features (5) and (6) are lo¬ 
cated closer in than the two arms and the submillimeter cavity, 
at previously inaccessible radii close to the star. Feature (5) may 
trace the optically thick innermost disk ( jlsella et aL||2008| , or 
alternatively, the inner arm of the SE spiral. 

Considering that the best-fit aspect ratio from the spiral mod¬ 
els is difficult to reconcile with the expected temperatures in this 
disk, and that on the other hand, the disk mass seems to be too 
low to account for gravitational instability, additional observa¬ 
tions coupled with advanced modeling are needed to understand 
which processes these non-axisymmetric features may trace. In 
particular, complementary observations on similar spatial scales 
with ALMA, which will trace different optical depths in the disk, 
is very promising. Future DPI observations with SPHERE and 
its visible-light instrument ZIMPOL (Roelfsema et al. 2010[ ) are 
also expected to place new constraints on the innermost struc¬ 
tures. 
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Appendix A: U# image and SNR map 

Figure |AT| presents the U$ images. Since the scattered light from 
a circumstellar disk is expected to be linearly polarized in the az¬ 
imuthal direction under the assumption of single scattering, U$ 
contains no signal, but only noise of the same magnitude as the 
noise in the Q$ image. For multiple scattering, the assumption 
of polarization only in azimuthal direction still holds approxi¬ 
mately, especially for close to face-on disks, and the deviations 
from azimuthal polarizations are expected to be well below the 
noise level ( Ave nha us et al.[20141 ). 

Figure [A72] provides qualitative SNR maps. They were ob¬ 
tained by first smoothing both Q$ and U# with a Gaussian kernel 
of 37 mas (close to the instrument PSF) and then dividing the Q$ 
by the standard deviation in the U$ images. The standard devi¬ 
ation in the U# image was calculated over a Gaussian-weighted 
area centered on the respective image point with a FWHM of 260 
mas. We note that these do not strictly provide signal-to-noise ra¬ 
tio estimates. Indeed, the noise reported in these maps depends 
on the area considered for the calculation of the variance from 
the U(f, images. Moreover, the first step smoothes out the readout 
noise and other small-scale (smaller than the smoothing kernel) 
noise (but not large-scale systematic noise). These maps are still 
very useful, as a feature detected in the Q$ image and appearing 
at a location of the map that has a high SNR value, should be 
considered as real. An example is feature 4, which lies in an area 
of high SNR and can thus be considered real, even though it does 
not seem apparent from the SNR maps alone. 

Appendix B: Probability plot 

Figure |B.1| shows the 2D probability distributions for each pair 
of free parameters. 
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Figure A.l. images (left, December 2014; right, March 2015). East is toward the left. Each pixel has been scaled with the square 
of its distance from the star, r 2 , to compensate for the r -2 dependency of the stellar illumination. The color scale is arbitrary. 
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Figure A.2. SNR maps with contours every 3 <x (left, December 2014; right, March 2015) 
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Figure B.l. Bayesian probabilities for the simultaneous fit of the two spirals. The top panel of each column provides the integrated 
ID probability for each parameter. The lower panels in each column are the 2D probability maps for each pair of parameters. 
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